Time‐Varying Latent Effect Model for Longitudinal Data with Informative Observation Times

Summary In analysis of longitudinal data, it is not uncommon that observation times of repeated measurements are subject‐specific and correlated with underlying longitudinal outcomes. Taking account of the dependence between observation times and longitudinal outcomes is critical under these situations to assure the validity of statistical inference. In this article, we propose a flexible joint model for longitudinal data analysis in the presence of informative observation times. In particular, the new procedure considers the shared random‐effect model and assumes a time‐varying coefficient for the latent variable, allowing a flexible way of modeling longitudinal outcomes while adjusting their association with observation times. Estimating equations are developed for parameter estimation. We show that the resulting estimators are consistent and asymptotically normal, with variance–covariance matrix that has a closed form and can be consistently estimated by the usual plug‐in method. One additional advantage of the procedure is that it provides a unified framework to test whether the effect of the latent variable is zero, constant, or time‐varying. Simulation studies show that the proposed approach is appropriate for practical use. An application to a bladder cancer data is also given to illustrate the methodology.     

Embryonic development in ballan wrasse Labrus bergylta

Eight primary embryonic developmental stages were assigned to eggs of ballan wrasse Labrus bergylta using key morphological features following standardized nomenclature: Ia, Ib, II, III, IV, V, VI and VI+, reared from single family clutches under comparable environmental conditions in Ireland and Norway. Development in L. bergylta is typical of demersal marine finfish species with a short egg stage. Hatching occurred c. 123 h post-fertilization (hpf) equivalent to 62·5 degree days at 12·2 ± 1·10° C (mean ±s.d.), after which the larvae swam intermittently near the surface of the water column.     

Five hepatopancreatic and one epidermal chitinases from a pandalid shrimp (Pandalopsis japonica): cloning and effects of eyestalk ablation on gene expression

Six cDNAs encoding chitinase proteins in Pandalopsis japonica were isolated by using polymerase chain reaction (PCR) cloning methods and bioinformatic analysis of expressed sequence tags (ESTs). The cDNAs, designated Pj-Cht1, 2, 3A, 3B, 3C, and 4, encoded proteins ranging from 388 to 607 amino acid residues in length (43.61-67.62 kDa) and displayed a common structural organization: an N-terminal catalytic domain, a Thr/Pro-rich linker region, and either 0 (Pj-Cht2, 3A), 1 (Pj-Cht1, 3B, and 3C), or 2 (Pj-Cht4) C-terminal chitin-binding domain(s) (CBD). Pj-Cht1 and 2 lacked the 5′ end of the open reading frame (ORF); the other Pj-Chts contained the complete ORF. All known decapod crustacean chitinases were segregated into at least four groups based on phylogenetic analysis and domain organization. Group 1 chitinases, represented by Pj-Cht1, were most closely related to insect group I chitinases and may function in the digestion of the peritrophic membrane. Group 2 chitinases including Pj-Cht2 show different domain organizations and pI value from other chitinases and appear to function in degradation of the old exoskeleton during the premolt period. Group 3 chitinases, represented by Pj-Cht3A, 3B, and 3C, may function in digestion of chitin-containing food and defense against pathogens. Group 4 chitinases, represented by Pj-Cht4, have two CBDs and their functions are unknown. Five Pj-Chts (Pj-Cht1, 3A, 3B, 3C, and 4) are expressed in the hepatopancreas and intestine, whereas Pj-Cht2 is expressed in epidermis and SG/XO complex suggesting crustacean chitinases can be classified into two groups (hepatopancreatic and epidermal) based on the expression profile. Eyestalk ablation (ESA) down-regulated the hepatopancreatic chitinase expression (Pj-Cht1, 3A, and 3C); Pj-Cht3B expression was not significantly affected by ESA. By contrast, mRNA levels of Pj-Cht2 were significantly upregulated in 7 days post-ESA. Pj-Cht4 mRNA levels were too low for measurement with quantitative polymerase chain reaction. ESA had no significant effect on chitinase expression in the intestine. These data indicate that Pj-Cht1, 3A, 3B, 3C, and 4 are hepatopancreatic chitinases that may function in the digestion of ingested chitin and the modification of peritrophic membrane in the intestine. By contrast, epidermal chitinase, Pj-Cht2 may play a role in chitin metabolism during molt cycle as shown in other crustacean group 2 chitinases.     

Around the world in 10 million years: biogeography of the nearly cosmopolitan true toads (Anura: Bufonidae)

Aim The species‐rich family of true toads (Anura: Bufonidae) has been the focus of several earlier studies investigating the biogeography of geographically widespread taxa. Herein, we employ newly developed Bayesian divergence estimate methods to investigate the biogeographical history of this group. Resulting age estimates are used to test several key temporal hypotheses including that the origin of the bufonid clade pre‐dates Gondwanan vicariance (~105 million years ago, Ma). Area cladograms are also invoked to investigate the geographical origin of the family. Location Worldwide, except the Australia–New Guinea plate, Madagascar and the Antarctic. Methods A phylogenetic hypothesis of the relationships among true toads was derived from analysis of 2521 bp of DNA data including fragments from three mitochondrial (12S, tRNA^val, 16S) and two nuclear (RAG‐1, CXCR‐4) genes. Analysis of multiple, unlinked loci with a Bayesian method for estimating divergence times allowed us to address the timing and biogeographical history of Bufonidae. Resulting divergence estimates permitted the investigation of alternative vicariance/dispersal scenarios that have been proposed for true toads. Results Our area cladogram resulting from phylogenetic analysis of DNA data supports a South American origin for Bufonidae. Divergence estimates indicate that the family originated earlier than had been suggested previously (78–99 Ma). The age of the enigmatic Caribbean clade was dated to the late Palaeocene–early Eocene. A return of bufonids to the New World in the Eocene was followed by rapid diversification and secondary expansion into South America by the early Oligocene (Rupelian). Main conclusions The South American origin of Bufonidae in the Upper Cretaceous was followed by relatively rapid expansion and radiation around the globe, ending with a return to the Americas via a Eurasian/North American land bridge in the Eocene. Though the exact route of this dispersal (Beringia or North Atlantic) remains unclear, an argument is made for the less frequently invoked North Atlantic connection. The origin of the enigmatic Caribbean lineage was found to be consistent with colonization following the bolide impact at the K/T boundary. These findings provide the first, firm foundation for understanding true toad divergence times and their truly remarkable and global radiation.     

Dielectric behavior and atomic structure of serum albumin

After 1946, serum albumin was available for studies. Its residue sequence and internal disulfide bonding was developed by 1976. We began to make dielectric dispersion studies and apply Perrin's equations for rotational relaxation times around the two axes of revolution in 1938. These data indicated that albumin should have an elongated shape. In 1992 atomic structure data indicated the molecule was heart-shaped. A similar 1998 study of albumin complexed with fatty acid showed that the molecule was substantially rearranged. We found that the dielectric constant of albumin solutions was sensitive to fatty acid content, making this property an attractive probe in stop–flow kinetic studies. Such studies show that the fatty acid reaction is a two-step process. The fatty acid first binds to exterior sites in a diffusion-limited second order reaction complete in 1 ms. Then a first order rearrangement reaction with 400 ms half-life follows. Thus the highly specialized serum albumin sequence of amino acid residues determines not only the structure of the unligated molecule, but also the distinctive structures of the numerous multiligated molecules.     

Comparison of (13)C(alpha)H and (15)NH backbone dynamics in protein GB1

This study presents a site-resolved experimental view of backbone C(alpha)H and NH internal motions in the 56-residue immunoglobulin-binding domain of streptococcal protein G, GB1. Using (13)C(alpha)H and (15)NH NMR relaxation data [T(1), T(2), and NOE] acquired at three resonance frequencies ((1)H frequencies of 500, 600, and 800 MHz), spectral density functions were calculated as F(omega) = 2omegaJ(omega) to provide a model-independent way to visualize and analyze internal motional correlation time distributions for backbone groups in GB1. Line broadening in F(omega) curves indicates the presence of nanosecond time scale internal motions (0.8 to 5 nsec) for all C(alpha)H and NH groups. Deconvolution of F(omega) curves effectively separates overall tumbling and internal motional correlation time distributions to yield more accurate order parameters than determined by using standard model free approaches. Compared to NH groups, C(alpha)H internal motions are more broadly distributed on the nanosecond time scale, and larger C(alpha)H order parameters are related to correlated bond rotations for C(alpha)H fluctuations. Motional parameters for NH groups are more structurally correlated, with NH order parameters, for example, being larger for residues in more structured regions of beta-sheet and helix and generally smaller for residues in the loop and turns. This is most likely related to the observation that NH order parameters are correlated to hydrogen bonding. This study contributes to the general understanding of protein dynamics and exemplifies an alternative and easier way to analyze NMR relaxation data.     

Genetic variation in organisms with sexual and asexual reproduction

The genetic variation in a partially asexual organism is investigated by two models suited for different time scales. Only selectively neutral variation is considered. Model 1 shows, by the use of a coalescence argument, that three sexually derived individuals per generation are sufficient to give a population the same pattern of allelic variation as found in fully sexually reproducing organisms. With less than one sexual event every third generation, the characteristic pattern expected for asexual organisms appear, with strong allelic divergence between the gene copies in individuals. At intermediary levels of sexuality, a complex situation reigns. The pair-wise allelic divergence under partial sexuality exceeds, however, always the corresponding value under full sexuality. These results apply to large populations with stable reproductive systems. In a more general framework, Model 2 shows that a small number of sexual individuals per generation is sufficient to make an apparently asexual population highly genotypically variable. The time scale in terms of generations needed to produce this effect is given by the population size and the inverse of the rate of sexuality.     

Identification of QTL for growth- and grain yield-related traits in rice across nine locations of Asia

Rice double-haploid (DH) lines of an indica and japonica cross were grown at nine different locations across four countries in Asia. Genotype-by-environment (G x E) interaction analysis for 11 growth- and grain yield-related traits in nine locations was estimated by AMMI analysis. Maximum G x E interaction was exhibited for fertility percentage number of spikelets and grain yield. Plant height was least affected by environment, and the AMMI model explained a total of 76.2% of the interaction effect. Mean environment was computed by averaging the nine environments and subsequently analyzed with other environments to map quantitative trait loci (QTL). QTL controlling the 11 traits were detected by interval analysis using mapmaker/qtl. A threshold LOD of >/=3.20 was used to identify significant QTL. A total of 126 QTL were identified for the 11 traits across nine locations. Thirty-four QTL common in more than one environment were identified on ten chromosomes. A maximum of 44 QTL were detected for panicle length, and the maximum number of common QTL were detected for days to heading detected. A single locus for plant height (RZ730-RG810) had QTL common in all ten environments, confirming AMMI results that QTL for plant height were affected the least by environment, indicating the stability of the trait. Two QTL were detected for grain yield and 19 for thousand-grain weight in all DH lines. The number of QTL per trait per location ranged from zero to four. Clustering of the QTL for different traits at the same marker intervals was observed for plant height, panicle number, panicle length and spikelet number suggesting that pleiotropism and or tight linkage of different traits could be the possible reason for the congruence of several QTL. The many QTL detected by the same marker interval across environments indicate that QTL for most traits are stable and not essentially affected by environmental factors.     

Heat capacities and a snapshot of the energy landscape in protein GB1 from the pre-denaturation temperature dependence of backbone NH nanosecond fluctuations

Protein stability is usually characterized calorimetrically by a melting temperature and related thermodynamic parameters. Despite its importance, the microscopic origin of the melting transition and the relationship between thermodynamic stability and dynamics remains a mystery. Here, NMR relaxation parameters were acquired for backbone 15NH groups of the 56 residue immunoglobulin-binding domain of streptococcal protein G over a pre-denaturation temperature range of 5-50 degrees C. Relaxation data were analyzed using three methods: the standard three-Lorentzian model free approach; the F(omega)=2omegaJ(omega) spectral density approach that yields motional correlation time distributions, and a new approach that determines frequency-dependent order parameters. Regardless of the method of analysis, the temperature dependence of internal motional correlation times and order parameters is essentially the same. Nanosecond time-scale internal motions are found for all NHs in the protein, and their temperature dependence yields activation energies ranging up to about 33kJ/mol residue. NH motional barrier heights are structurally correlated, with the largest energy barriers being found for residues in the most "rigid" segments of the fold: beta-strands 1 and 4 and the alpha-helix. Trends in this landscape also parallel the free energy of folding-unfolding derived from hydrogen-deuterium (H-D) exchange measurements, indicating that the energetics for internal motions occurring on the nanosecond time-scale mirror those occurring on the much slower time-scale of H-D exchange. Residual heat capacities, derived from the temperature dependence of order parameters, range from near zero to near 100J/mol K residue and correlate with this energy landscape. These results provide a unique picture of this protein's energy landscape and a relationship between thermodynamic stability and dynamics that suggests thermosensitive regions in the fold that could initiate the melting process.     

Autotetraploid tangor plant regeneration from in vitro Citrus somatic embryogenic callus treated with colchicine

In vitro chromosome doubling of embryogenic callus lines of the Citrus cultivars Umatilla and Dweet tangors (Citrus reticulata Blanco×C. sinensis [L.] Osb.), Caffin clementine (C. clementina Hort. ex Tan.) and Wheeny grapefruit (C. paradisi Macf) was carried out in the presence of either 0.05 or 0.1% colchicine, or 0.01, 0.05 or 0.1% oryzalin. Embryogenic callus development was partly suppressed in the presence of colchicine, and completely suppressed by oryzalin at all concentrations tested. No plants were regenerated from any of the oryzalin treatments. Ploidy level of plants regenerated from the colchicine treatments was determined using flow cytometry and chromosome squashes. Three desirable non-chimeric, autotetraploid plants of the mono-embryonic cultivar Umatilla were produced using 0.05% colchicine and one from 0.1% colchicine. One mixoploid Dweet plant was produced using 0.1% colchicine.